Conserved superenergy currents 



Ruth Lazkoz 1 *, Jose M.M. Senovilla 1§ , and Raiil Vera 21 
1 Ffsica Teorica, Universidad del Pais Vasco, 
Apartado 644, 48080 Bilbao, Spain 
2 School of Mathematical Sciences, Queen Mary, University of London, 
Mile End Road, London El 4NS, U.K. 

February 7, 2008 



Abstract 

We exploit once again the analogy between the energy-momentum tensor and the so- 
called "superenergy" tensors in order to build conserved currents in the presence of Killing 
vectors. First of all, we derive the divergence-free property of the gravitational superenergy 
currents under very general circumstances, even if the superenergy tensor is not divergence- 
free itself. The associated conserved quantities are explicitly computed for the Reissner- 
Nordstrom and Schwarzschild solutions. The remaining cases, when the above currents 
are not conserved, lead to the possibility of an interchange of some superenergy quantities 
between the gravitational and other physical fields in such a manner that the total, mixed, 
current may be conserved. Actually, this possibility has been recently proved to hold for 
the Einstein-Klein-Gordon system of field equations. By using an adequate family of known 
exact solutions, we present explicit and completely non-obvious examples of such mixed 
conserved currents. 



1 Introduction and basic results 



Conserved quantities in Gravitation are either defined only at "infinity", see e.g. and 
references therein, or they arise as a consequence of the existence of divergence-free vector fields, 
called local currents, and of the use of Gauss' theorem applied to appropriate domains of the 
spacetime, see e.g. [1'dl 1221 l2"rjj . A traditional way of constructing divergence- free vector fields 
is by using the intrinsic symmetries of the spacetime 1 : if £ is a Killing vector field, V( a £j3) = 
0, representing an infinitesimal isometry, then many local currents and associated conserved 
quantities arise, such as for instance every vector field of the form A£, where A is any function 
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X A spacetime (V, g) is a 4-dimensional manifold V with a metric g of Lorentzian signature (-,+,+,+). Indices 
in V run from to 3 and are denoted by Greek small letters. Square and round brackets denote the usual (anti-) 
symmetrisation of indices. The tensor and exterior products are denoted by ® and A, respectively, boldface 
letters are used for 1-forms and arrowed symbols for vectors, and the exterior differential is denoted by d. The 
Lie derivative with respect to the vector field £ is written as and the covariant derivative by V. Equalities by 
definition are denoted by =, and the end of a proof is signalled by M. 
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satisfying £(A) = 0. This is a simple consequence of the invariance of the geometrical background 
under the Killing vector field. Of course, if there are several Killing vectors £l m the spacetime, 
then any linear combination of the vector fields of type A l £l is also divergence-free whenever 
each A L is invariant with respect to the corresponding 

There are also other, physically more interesting, types of locally conserved currents con- 
structible from £ together with other objects. Take, for example, the typical case of for a 
divergence- free symmetric tensor t^ u . Obviously we have 

v M (#•"&,) = (V) s„ + tw v M e„ = o 

where the Killing property V^^) = has been used. Observe that both the divergence-free 
property and the symmetry of t^ v are needed here. In fact the divergencelessness of t^ v is a 
sufficient but not necessary condition, as one only needs that £ and V-i be orthogonal for the first 
term to vanish. Furthermore, observe that just a conformal Killing vector field (V^^) oc g a p) 
is required if t^ u is also traceless. 

A particular important example of this type of locally conserved currents is provided by the 
case of the total energy-momentum tensor T^ u of the matter fields in General Relativity. In 
this case, due to the Einstein field equations, T^ u is symmetric and diver genceless, and thus any 
Killing vector field £ provides a divergence-free current defined by 

-ia/3 



J a (£) = T^Cp. (1) 

These currents have also a physical meaning which depends on the particular character of the 
Killing vector £. For instance, in stationary spacetimes there is a timelike Killing vector £, and 
the corresponding current J (^j can be thought of as the local energy-momentum vector of the 
matter in the stationary system of reference. Similarly, for translational or rotational Killing 
vectors one obtains currents representing the whole linear or angular momentum, respectively, 
at least in flat spacetimes. 

Two important remarks are in order here. First, in many cases the mentioned currents J (^j 

are of the form referred to above, that is, J ( £ ) = B l £l for some particular functions B L . This 



does not deprive J (^j of its physical significance. Rather one has to see this result as stating 

that, among the huge variety of divergence-free vectors of the form A l £l with general A L , the 
particular one with the functions A L = B L have a specially relevant physical meaning. Several 
cases in which this happens are given by the following known result [T%| 12*3 } FH fTT>t ITTT] . 2 

Result 1.1 Assume Einstein's field equations of General Relativity hold. 
1. If £ is a hyper surf ace- orthogonal Killing vector, then *7" (£) A £ = 0. 



2. If ^ and ff are two commuting Killing vectors [£,ff\ = and they act orthogonally transi- 
tively on non-null surfaces, then J (^j A £ A T] = 3 (ff) A £ A rj = 0. 

2 Recall that a vector field v is called hypersurface-orthogonal (also "integrable" ) if v A dv — 0, and that two 
vector fields orthogonal to the non-null surfaces spanned by two given vector fields v and w generate surfaces 
whenever the two 1-forms v and w satisfy v Aw A dw = v A w A dv = 0. Thus a G2-group of motions generated 
by two Killing vectors £ and ff is said to be acting orthogonally transitively on non-null orbits when there exists 
a family of surfaces orthogonal to the orbits, which means that £A77Ad?7 = £A77Ad£ = 0. 
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The first of these results means that if there are (locally) hypersurfaces orthogonal to the Killing 
vector, then the corresponding current lies in the direction of the Killing vector itself. This 

implies that the integration of J {t\ •£ (and of J f|j over any two different hypersurfaces 

orthogonal to £ is the same. The second one states that any two (group-forming) Killing vectors 
whose surfaces of transitivity are orthogonal to a family of surfaces have currents which lie 
necessarily in the orbits of the G2 group generated by them. Of course, both results are quite 
natural and follow from the invertibility of the orthogonally transitive Abelian G2 group (which 
reduces to hypersurface-orthogonality for 1-dim groups) 

The second remark concerns the case when there are several types of matter fields in the 
spacetime, and the corresponding energy-momentum tensors. For simplicity, assume that there 
are two different sorts of matter contents (an electromagnetic field and a fluid for instance), with 
corresponding energy-momentum tensors T^} and T^), forming an isolated system. Then, only 

the total energy-momentum tensor T a p = T^l +T^) is divergence- free, and therefore the partial 

currents constructed from T^J or T£a are not conserved separately in general. Of course, this is 
physically meaningful and reasonable, as it leads to the conservation of the corresponding total 
quantity, and thereby implies the exchange of some energy-momentum properties between the 
two fields involved. A classical example of this situation can be found in the standard reference 
19 for the case of Special Relativity, where one can see that if the sources of the electromagnetic 
field are taken into account then the energy-momentum tensor of the electromagnetic field is not 
divergence-free, but its combination with the energy-momentum tensor of the charges creating 
the field provides a total energy-momentum tensor which is divergence-free and thereby provides 
the usual conserved currents associated to the ten Killing vectors of flat spacetime. This proves in 
particular that the energy-momentum properties can be transferred from one field to another, 
or from matter to the fields, and vice versa. This transference property together with the 
positivity and the conservation are the basis for the paramount importance of the concept of 
energy-momentum in physical theories. 

In this paper we shall analyse other locally conserved currents arising in Gravitation which 
keep the mentioned three properties: positivity, local conservation, and interchange. These local 
currents arise due to the existence of the so-called "superenergy" tensors, see PJEHHHUEHIEZIEE] 
and references therein, and, as we shall prove, they have analogous attributes to those so far 
mentioned for the energy-momentum currents (jljl. The type of currents we shall analyse have 
expressions such as 

or adequate generalisations of this using different Killing vectors, where T a ^ Xfl will be any of 
the available superenergy tensors for the physical fields involved. 

In particular, we start by constructing the corresponding currents for the Bel superenergy 
(s-e) tensor (2J E3 HE] > which is divergence- free in vacuum (where it is called the Bel- Robinson 
tensor) and in Einstein spaces, and, therefore, it can lead to conserved s-e currents. We shall 
see in Section |21 that, in fact, all the currents constructed from the Bel tensor (in vacuum or 
not) and a Killing vector satisfy properties analogous to those of Result II. 1| which will result 
in divergence- free currents in many general situations. Some specific examples of these currents 
are written down and analysed in Section |3] leading to physical consequences of some interest. 

Nevertheless, in completely general situations, the Bel s-e currents will not be divergence- 
free; something which seems natural if the s-e concept is to have any physical meaning at all. 
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This is because if in the spacetime other fields than the gravitational one are present, one should 
expect in principle that the conserved current be formed by a combination of the sum of the 
gravitational Bel s-e current and the s-e current of the matter fields involved, in analogy with 
the energy-momentum currents, as explained before. This is supported also by the fact that 
the s-e currents of the matter fields are conserved in flat spacetime -that is, in the absence of 
gravitational field-, see for a proof |29l 1281 15U] , As a matter of fact, the existence of such mixed 
divergence- free currents has been rigorously proved in |29|l28 for the case of a minimally coupled 
scalar field, that is, for the Einstein-Klein-Gordon situation, whenever there is a Killing vector. 
As explicitly shown in |291 128j . none of the two single currents, the scalar field s-e current nor 
the Bel one, are divergence-free separately in general. These results are summarised briefly in 
Section 0J 

One of the purposes of this paper is to compute these divergence-free mixed currents in 
some particular situations, thereby providing explicit examples of the exchange of s-e between 
the scalar and gravitational fields, and the conservation of combined quantities. However, there 
arise several technical difficulties. First of all, the result for the Bel tensor which is analogous 
to Result 11.11 above, presented in Results 12.11 and 12.21 of Section implies that (see Theorem 
12.1(1 : 1) if the Killing vector is hypersurface-orthogonal, then its corresponding Bel current 
is automatically divergence- free; and 2) if two Killing vectors form an Abelian group which 
acts orthogonally transitively on non-null surfaces, then all the corresponding Bel currents are 
divergence-free too. Thus, even though these results are physically interesting on their own, 
they prevent the existence of mixed currents in those cases. Hence, we have to resort to other 
more general cases in order to find such mixed currents. But then a second problem arises, 
namely, to find explicit spacetimes which are solution of the Einstein-Klein-Gordon equations 
under the needed more general circumstances: there are hardly any known solutions with a 
non-hypersurface-orthogonal Killing vector, and very few with a non-orthogonally transitive 
2-parameter group. 

The good news is, though, that there is at least an explicitly known family of solutions of 
the type we need. These solutions were found by Wils in |34| and we will use them in the second 
part of Section 0] to construct the sought mixed divergence-free s-e currents. This will provide 
explicit examples of the s-e interchange and of conserved currents which would have been very 
difficult to find, and even more to be singled out, if the concept of superenergy tensors had not 
been available. 

2 General results on Bel currents 

We will use the notation and terminology of 28 for the superenergy tensors. The Bel tensor, 
introduced by Bel more than 40 years ago |H is the basic superenergy tensor of the total 
gravitational field, that is to say, the s-e tensor T a p\^ {-R[2],[2l} associated to the Riemann tensor 
Ra/3Xp- Rs explicit expression reads [TU1 151 125] 

BafiXfi = RapXoRjf p° + Rapno-R(3 P \ a ~ 
1 1 1 

— ydapRprXo-R^ p" — -^9XtJ,RapaTRf3 PaT + -g a /39XpRpTCTuR pTaL ' , (2) 

from where the following properties explicitly arise: 

B a f3Xp, = B(a/3)(Xp) = B X p a f3, B P pXlJl = 0. (3) 



4 



Using the second Bianchi identity V \ v R a p]x^ = one obtains the following expression for 
the divergence of the Bel tensor: 

V Q B a ^ = R% x a J^ p + RP/ a J Xa P - ^R^J^p, (4) 

where Jxu/3 = —J^x/3 = — and Rp^ is the Ricci tensor. Note that because of (0) 

this is the only independent divergence of the Bel tensor. From Q) we obtain the fundamental 
result that B is divergence-free when Jxap = 0. This happens, for instance, in all Einstein spaces 
(where R^ u = Ag uu ), so that in particular this implies that the Bel- Robinson tensor [Tl 01 HI 05]. 
which is the s-e tensor of the Weyl conformal tensor and coincides with B in Ricci-flat spacetimes 
jS], is divergence- free. The divergence-free property of the Bel or Bel- Robinson tensors in these 
cases allows for the construction of divergence-free currents whenever there is a Killing vector. 

These currents are built in a similar way to those formed with the energy-momentum tensor, 
as we show next. Following |2H] one can define the Bel current with respect to any three Killing 
vector fields £ 2 > £3 (they do not need to be different!) as 



fa (ll> ^2; £3) = B (aP\)n fl £2 ^3 = B (a/3X f j,) fl ^2^3- ( 5 ) 

To avoid unnecessary writing we will omit the repetition of the Killing vectors: j (^,1,6,2^ = 



3 [£i,£,2,£,2] and j (£) = j The divergence of these currents can be easily computed 

to give 

where we have used the fact that the £ a are Killing vectors: V (a ^f } = (a = 1,2,3.) Therefore, 

the vanishing of the divergence of the Bel tensor implies the vanishing of that of j \ £2, £3) , 

defining, thereby, a conserved current. 3 

Nevertheless, the Bel currents are not only conserved when the Bel tensor is divergence- 
free. Following the general arguments in [Hj showing the invertibility of orthogonal transitive 
Abelian G2 groups (and, in particular, integrable 1-dim groups) acting on non-null orbits, since 
the Bel tensor is intrinsically defined, if the metric admits an invertible group then the Bel 
tensor must be invertible implying that its contraction with any odd number of Killing vectors 
provides a tensor tangent to the group orbits. Therefore, the s-e currents constructed from a 
hypersurface-orthogonal Killing vector or two non-null Killing vectors generating an Abelian 
orthogonally-transitive G2 will necessarilly be tangent to the orbits of the corresponding groups. 
This is stated as follows: 

Result 2.1 If £ is a hypersurface-orthogonal Killing vector, then j A £ = 0. 

Result 2.2 Let £ and ff be two independent commuting Killing vector fields spanning non-null 
surfaces. Then, if they act orthogonally transitively, their four associated Bel currents satisfy 

j (t\A^Arj = j(fJ)A^Ai 1 =j (i,ff)A^Ar t =j (ff,£) A £ Arj = 0. (6) 



3 Observe that only in strictly Ricci-flat spacetimes we can use conformal Killing vectors for £ 1; £ 2 > £3 obtaining 
also local currents of type JKJ , as only in that case the Bel-Robinson tensor is divergence- and trace- free. Then, 
one can involve conformal Killing tensors — or even conformal Yano-Killing tensors as in |15|-. 
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These currents satisfy then properties completely analogous to that of the energy-momentum 
currents stated in Result IT~T1 and therefore they will be divergence-free in quite general cases, 
as we will show below. For the sake of completeness, we provide the identities involving the s-e 
currents and the geometrical properties of 1-dim and 2-dim groups of isometries, from which 
Results O and O trivially follow. 

Identity 2.1 Let £ be a Killing vector field in (V,g). Then we have 

e^^BapXl^y] = I {Cx^R^Vr (£ [M £ 7;CT] ) + ^R^^Va (^ ;r] ) } . (7) 

Proof. Since £ is a Killing vector we have j!21 122j 

V«V/9& = ( p R P a(SX, (8) 

so that direct contraction of the Bel tensor © with £ yields, on using (JHJ), 
The exterior product of this one-form with £ gives then 

ef^B aPX ^ y] = 26vv ff e A v p v (T e [ ^ 7 ] + ^av^v^^^. (9) 

Expanding the covariant derivative of the three-form £ A d£ we get 

3 V P (£[ 7 6x; M ]) = -3 V 7 £ [p V^ ff] + kVpV^ - 2 VpV ff ^er] • ( 10 ) 

and taking into account that CaV^V 7 ^ is symmetric due to (JHJ), its contraction with (|1UI) 
provides 

3 eAv^v CT e A v p (e [7 e CT;M ]) = -2 6vvYv p v^ 7] . (n) 

Let M^> pp,a be any tensor with an arbitrary number of indices (denoted by {£1}) plus three 
indices such that M^^ pfJj<7 ^ = M^ p ^ a ^ = 0. After relabelling the indices, the contraction of 
(fm|l with AfWw and £ 7 gives 

3 e /? M {^} CT pr V(T (£ [/3 £ T;p] ) = e A ^V CT V p V r M^>^ + 2 ^V ff V T ?^ p M{"^^. (12) 

Replacing now using (JHJ and (jll|) we find 

3 ^^Vt]^ (e[^r; P ]) = V ff V P £ T R P ™ ^ - 3 &V'V ff £ A V p fe^ iA4 ]) . (13) 
The final step is just to compute the linear combination (|9*|)+ Q11[1 — h (J13|) to obtain 

which, after using (JHJ) and some index manipulation, gives the desired identity (J7J). ■ 

Note that i a eeB a pxp = C^ X B^ x)lt = t a ^ t X B [a(3Xli) due to the symmetries in ©. If 
£ is hypersurface-orthogonal, i.e. (,[ p (,f3 ;fl ] = 0, Result l2~T1 explicitly follows at once as a corollary 
of formula (|7j). 
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Identity 2.2 Let £ and fj be two Killing vector fields in (V, g) generating a two dimensional 
group of isometries acting on non-null surfaces. Taking into account i3ty . \4$ , \3b )) -\3y \) in 

the Appendix, the four associated currents can be expressed as follows 

-3W 2 j [a (£a,£b) «V] = 

+ {^U^^BfiW 60 + ^^t^A/ju^) ("aSi p - i x t 2p )R T P t[a w^] + 

+ \ [(2 (6l • U)£bt + 4(^4 • £b)IUt) ^ + 



l (2w X(7 + *w X(7 ) w^R^xa ((U ■ U)& + 2(6i • (%S 1[a u; H - ^S 2[Q u; H ), (14) 



where the indices A,B,C take the values 1,2 to denote £ and 77 respectively, and * denotes the 
usual Hodge dual. 

Proof. Not to overwhelm the text, we prefer to present the proof in an Appendix. ■ 

Clearly, if the G2 isometry group is Abelian and acts orthogonally transitively, then S^Atu = 
Q,A/\w = and ~Ea = = for A = 1, 2, as follows from ((36|) - (|39|) . and Result l2~2l is recovered. 

Let us introduce the notation jy, (T = 1, 2, 3, 4), for the four Bel currents appearing in ©. 
Then, Result l2~2*l implies that 

jr = a T {x)i + b r (x)ff. (15) 

The functions ax and by are not arbitrary. Note that, given that £ and fj are Killing vectors, 
then £^R a px P = £f]R a p\n = 0, and, therefore, 

£gB a p\n = £tfB a px\x = 0. (16) 

Hence, if we compute the Lie derivative with respect to both Killing vector fields of equation 
(|T5|) , by using (fl6|) and the definition © of Bel currents we deduce 

= £fr = (£g a T ) £+ (£g by) fj 

and a similar relation replacing £^by £fj. In consequence, the functions ay and by are restricted 
by 

£ 1* ay = £g by = £fj ay = £^ by = 0. 

Taking now the divergence of equation Q15|) and using this result we arrive at 

V p j T = T = 1,2,3,4. 

The case with only one hypersurface-orthogonal Killing vector can be treated as a special case 
of equation 1)15(1 with T = 1 and b\ = 0. All these results can be summarised as follows: 

Theorem 2.1 1. If £ is a hypersurface-orthogonal Killing vector then its corresponding Bel 
current j ( £ ) is proportional to £ and divergence- free. 
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2. If £ and fj are two commuting Killing vectors acting orthogonally transitively on non-null 
surfaces, then their four corresponding Bel currents jx He in the 2-planes generated by 
{£, fj} and are divergence-free. 

These results are very interesting, and reinforce, once more, the various analogies between the 
energy-momentum tensors and some of the superenergy tensors. Nevertheless, as remarked 
in the Introduction, they forbid the existence of mixed superenergy currents if they are to be 
constructed with Killing fields of the type appearing in the results. Unfortunately, most of the 
known Einstein-Klein-Gordon exact solutions contain such types of Killing vector fields. In order 
to get a flavour of the possible interpretation of the s-e currents, we are going to present some 
simple cases of pure divergence-free currents in the next section |3J leaving the more difficult 
mixed case for section |IJ in which an explicit solution of the field equations without the above 
mentioned properties for the Killing vector fields will be used. 

3 Bel and Bel-Robinson currents in Reissner-Nordstrom and 
Schwarzschild spacetimes 

Let us consider the Reissner-Nordstrom spacetime, whose line-element, in standard spherical 
coordinates {i, r, 9, ip}, reads [TU1 l2*2l ITH] 



ds 2 = - ( 1 - — + 41 ^ 2 + f 1 - — + 4 ) - '■W- ( l7 » 



where m and q are arbitrary constants representing the total mass and electric charge of the 
particle creating the gravitational field, and d£l 2 = d6 2 + sin 2 6 dip 2 is the canonical line-element 
in the unit 2-sphere. We only consider the exterior asymptotically flat region defined by 1 — 
+ ^ > which restricts the range of the coordinate r to r > if m 2 < q 2 , or to r > r + = 
m + a/ m 2 — q 2 if m 2 > q 2 . The metric (|17j) is the general spherically symmetric solution of the 
Einstein-Maxwell equations. The electromagnetic field is given by F = (q/r 2 ) dt A dr, and its 
energy-momentum tensor expressed (in natural units 8irG = c = 1 where c is the speed of light 
in vacuum and G is the gravitational constant) in the coordinate basis is 

q 2 

T^ = ^xdiag {-1,-1, 1,1}. (18) 

The particular case with q = is the vacuum Schwarzschild solution (then r > 2m). 

In the chosen range of coordinates, the metric is static (£ = d/dt is a hypersurface-orthogonal 
timelike Killing vector) and spherically symmetric. We enumerate the four Killing vectors £l, 
(L = 1,2,3,4) as 

& = W & = ™<P-jjg+<**<P«*0^, f 3 = cos^-sin^cot^, U = ^. (19) 

We want to study the Bel (and Bel-Robinson for the q = 0-case) currents constructed with 
the Killing vectors (|19|) for these spacetimes, together with the currents arising from the energy- 
momentum tensor. These will lead to some conserved quantities. Of course, one expects that 
all the conserved quantities will be just functions of the constants m and q, but our aim is to 
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ascertain which particular combinations of m and q arise at the energy and superenergy levels. 
To find them we compute explicitly all the Bel currents: 



. • 1 / 2m q 

m) = -- ( i — + 



Q 2 \ 2 <f 
6 ( m + -jr 



(20) 



. 1 / 2m q 2 



2\ 2 



12 ( m - — + \ 



3 (faj&jfc) = ^6 



12 ( m - — ^ 



(21) 



,2\ 2 



6 ( m + 



{ (ia ■ Cb) £c + (6 • fc) + (& • fa) , 



where a,b,c G {2, 3, 4}. Similarly, the energy-momentum currents defined in (J^) can be computed 
on this metric and they read: 



j(ia 



(22) 



It can be easily checked that all the above currents are divergence- free, and therefore they lead 
to a set of conserved quantities via Gauss' theorem applied to appropriate compact 4- volumes of 
the spacetime. For illustration purposes, we are going to compute now these conserved quantities 
for the simple but relevant case in which the compact region is taken to be bounded by two 
t =const. hypersurfaces, say t = t\ and t = t 2 , t 2 > t\, and two r = const, hyper surf aces, say 
t = T\ and r = r 2 , r 2 > r\. We call this region K. Applying Gauss' theorem we immediately get 



0. 



where A represents any of the j or J computed above, dK denotes the boundary of K, n 
is the outward unit normal to K and d 3 <7 is the canonical volume 3-element on dK. Clearly, 
dK = {t = ti, n<r<r 2 }U{t = t 2 , n < r < r 2 } U {r = n, h < t < t 2 } U {r = r 2 ,h < t < t 2 }, 
so that the corresponding unit normals are proportional to dt for the first two regions, and to 
dr for the remaining two. From the explicit expressions of the currents, we know that none of 
them has a non-zero component along dr, and therefore the integrals on {r = r{\ U {r = r 2 } 
vanish. Thus, only the integrals on {t = t±} and {t = t 2 } remain, and thus they must be equal. 
In other words, the integrals 

/ A"n M d 3 c7, (23) 

Jt=const. 

taken over any portion of a t =const. hypersurface bounded by the values r\ and r 2 , are constant 
in the sense that they are independent of the particular t =const. hypersurface. It is quite 
remarkable that the integrand in (|23j) for the case of the electromagnetic current IJ22JI is not the 
standard expression for the electromagnetic energy density, taken usually as T^ v n^n u . However, 
(|23l) leads to simpler results and to expressions which are apparently correct (at least in a naive 
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comparison with the analogous ones in the classical theory) such as ([25)) and (|26|) found later. 
Thus, one wonders if, in a stationary frame of reference defined by £, the correct expression for 
the electromagnetic energy density taken over an extended region of the spacetime — and not 
only at a point — should be in fact T^^^riy. 

Coming back to (|2"3)l . and by noting that n = — (l — "^r + % 

/ 2\-i/2 

I 1 ^ + 4 r 2 sin OdrdOdip, those integrals reduce to 



2 7r /"7T rr2 

dip / sin 6»cZ6» / (-A*)r 2 dr, (24) 

Jri 

where A 1 denotes the t-component of A. Observe that, once again, among the explicit expressions 
for j and J ', only those involving £1 have a non-zero t-component, and therefore they are the 
only ones which may give non-trivial conserved quantities: all other currents provide constants 
which are simply identically zero. Actually, not even all of the currents involving £1 provide a 
non-zero constant, because the angular integrals can also vanish. In summary, it is easy to check 
that the only non-trivial constants arise from the currents (|22|) , (|20j) and the three with a = b in 
(f2"T|) . these last three being, in fact, equal. Denoting them by Q±, Q\ and Q2 respectively, they 
clearly depend on the values of r% and r 2 . Their explicit expressions are: 

Ql =A^q 2 (---), (25) 



4 7T 

Qi = = (45m 3 r 3 - 90m 2 q 2 r 2 - 30m 2 r 4 + Qbmq A r + A5mq 2 r 3 - 21g 4 r 2 - 15g 6 ) 



15 r 7 

>■■: 

27r 3 



''2 



'1 



Q 2 = - (36m 2 r 2 - 36mq 2 r + 13q 4 ) . 

n 

It is also remarkable that these three expressions are strictly positive. This is a general result — 
independent of the spacetime if it is stationary — for the currents of type Q\ and Q%, and follows 
from the dominant energy condition satisfied by T^ v in the first case, and from the analogous 
dominant property satisfied by the Bel tensor j^U El El HE1 m the second, as the corresponding 
integrands are positive. This was noted for the Bel-Robinson case for instance in jllj (see also 
the recent generalisation in Nevertheless, there is no similar reasoning for Q2 and its 

positivity can only be inferred from the explicit expression of the local current (j21|) . 

If we want to obtain constants associated to the spacetime, we can also take the limit cases 
in which the above expressions are computed over a whole slice t =const. by taking r 2 — > 00 
and the minimum possible value for t\. There appear then two different situations depending 
on whether m 2 > q 2 or not. If m 2 < q 2 , then the minimum value for r\ is r\ =0, and clearly all 
the above constants diverge. This is natural and analogous to what happens in flat spacetime if 
we integrate the energy over the whole space. However, for the case with m? > q 2 , the existence 
of the event horizon at r + = m + y m 2 — q 2 provides a finite minimum value for n, and this 
leads to finite conserved constants. A simple computation leads to 



Q 



47r q 2 
r + 

Ait ( ,m , m 2 m 3 
6 - 14 h 14-5- - 5^- 

15 r+ V r+ ri. ri 
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Q 2 = —r+ 13-16— + 16 „ 
27 \ r + r+ 

For the particular case of the Schwarzschild spacetime (q = 0), they become simply 

„ „ 7r „ 16-7rm 

Si = 0, Qi = — , Q 2 = —3— • 

An obvious question arises: are these particular combinations of m and q special in any sense? 
And if yes, why? Clearly, the meaning of Qi is the total electrostatic energy of the Reissner- 
Nordstrom black hole with respect to the static observer. However, in order to check this simple 
statement one has to put back all the physical constants which were implicitly taken to be unity. 
Denoting by M and Q the total mass and charge in the correct units, respectively, we have 
m = GM/c 2 and q 2 = GQ 2 /c 4 . Now, the energy- momentum tensor (|18[) must be multiplied by 
c 4 /(87rG), and therefore we finally get for the physical Q\\ 

1 G 2 

Qi = I — ■ (26) 

2 r + 

This is a satisfactory result which may provide a physical interpretation for Qi, as this formula 
is reminiscent of that for the energy of the electric field in classical physics. 

Now, the same type of reasoning is needed for the superenergy quantities Q\ and Q 2 - To 
that end, we need to know which type of physical dimensions are carried by the superenergy 
tensors. This can be deduced from several independent works ^JEHSOlEl) and as is explained 
in [2Sj (p. 2820) the correct physical units for the physical Bel tensor seem to be energy density 
per unit surface, i.e. ML~ S T~ 2 . This means that we have to multiply the Bel tensor by kc^/G 
to get the correct physical quantities, where k is a pure number to be chosen. Doing that we 
can obtain the following expression for the physical Q\: 



Ql = K- T [ 7 Mc 2 + (3— )= K-j- t Mc 2 + 20Qi (27) 



4t/1,, 2 n Q 2 \ 4tt/1 

5-2- -M C 2 + /3— ]=K- T [- 

1/lm 2 23 m 4\ , ^ r __ 

where p = — — ^ 1 — is a dimensionless quantity. Of course, since (1271) is written 

2 \3 r+ 30 r + 5/ 

in terms of three non-independent quantities, M, Q and r + , it is not unique. It only provides 
a possible good choice. For the particular case of Schwarzschild's solution (q = =>■ r + = 2m) 
the previous quantity reduces simply to 

Qi = k-^Mc 2 . 

We note in passing that if we chose k = (2tt)~ 2 we could rewrite this as 

Qi = ^ (28) 

where A is the area of the horizon. Nevertheless, such an amusing or inspiring interpretation 
for Q\ is not possible in the general case (|2*7jl due to the factor (3. 

A possible way to try and obtain alternative, more complete, expressions of type (f2*7|) is by 
using the existence of independent superenergy of pure electromagnetic origin, see [2Sj. In order 



11 



to check this possibility, we can use any of the s-e tensors for the electromagnetic field (the most 
general one depends on six arbitrary constants, see |271 128j ) because the corresponding current 
constructed with is always the same, as can be easily checked. Thus, we can use the basic s-e 
tensor T a0Xfl {V [1} F [2] } for VF, given by [23I2E] 



E a p\n = VaFxpVpF/ + V a F w V p F x P - g a pV a F Xp X7 a F/ - 

—gxtNaFapVpF*" + \g a pgx l NrF (7p V T F' J P . (29) 

Its symmetry properties are 

E a (3\ti = #(a/3)(A,u) 

but it is not symmetric in the exchange of af3 with A// (a tensor with that symmetry was 
found many years ago by Chevreton [TQ1 EDj; Chevret on's tensor is simply proportional to 
F a /3X^i + E\u a p. In what follows, one can use either of these tensors or, for that matter, its 
completely symmetric part). The current associated to (|29]l. and analogous to (|20jl. can be 
easily computed to produce 

so that it is clearly divergence-free too. Thus, by a similar procedure as above it gives rise to a 
conserved quantity when integrated over any t =constant spatial section. Its value is 



jf) 47ro 2 / ,m 2 m \ 4tt ( m 3 m 2 _m 

})*> = — \ A— - 11 h 8 = 8-5- - 26^- + 27— 

35 ri \ rl r + / 35 r + V ri r± r+ 



Therefore, by forming linear combinations with positive coefficients c\ and ci of the respective 
s-e tensors we can get a conserved 
units have been restored, becomes 



s-e tensors we can get a conserved quantity given by c\Q\ + c 2 Q\ which, when the physical 



cxQx + c 2 Q[ F) = k% (aMc 2 + B— 



with 

r 2 

,,2 



Tfi TTl 

A = a + (b- 2a)— + {c x - 26) -5-, 



8 2 . 2 11 .m 2 4 ^m 2 

B = — c 2 - -ci + (a + —a - — c 2 )— + (6 - -ci + — c 2 — . 
35 5 15 35 r + 3 35 r+ 

Here a and b are spureous constants (due to the fact that 1 — 2m/r + + q 2 jrh =0), so that these 
expressions can be simplified on using this freedom. One can try to restrict the values of the 
constants by using other kind of arguments. For instance, there are independent results that 
prove the conservation of the simply added superenergy along shock-wave discontinuities for 
Einstein-Maxwell spacetimes |2U) (see also subsection 7.3 in [!2H] and references therein), hence 
it seems reasonable to put in any case c\ = c 2 . Setting then c\ = c 2 = 1 and, as an example, 
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choosing a and b such that A does not depend on m/r + so that the case Q = is explicitly 
recovered, we get 

r± V 4 



1 /ll m 2 29 m \ 

where now B = ( = h 6 ) . Other choices for a and b lead to similar formulas. 

35 V 6 r\ 12 r + J 

Thus, the interpretation of Q\ and, in general, of its combinations with Q\ , remains somewhat 
obscure. 

For the sake of completeness, we include now the physical Q2 (note that the corresponding 

( F) 

Q 2 vanishes, so that in this case the linear combinations of the superenergy tensors always 
lead to a quantity proportional to Q2): 

Q 2 = 8ttk (^Mc 2 + 7^ = 8vtk (5Mc 2 + 2 7 Qi) , 

where now 5 = 2(8/27 — k)m/r + + k + 10/27 and 7 = km/r + — 13/27 and k is a disposable 
constant similar to the previous a and b. Choosing again 5 not to depend on m/r+, one gets 



16vr 

Q2 = K — 



Mc + 9 ( r 



The units of this quantity are those of energy, that is ML 2 T~ 2 , but given its origin and con- 
struction one has doubts about whether considering it as such or as a kind of "supermomentum" 
per unit surface: (M L 2 T~ 2 )L 2 / L 2 . In the Schwarzschild case, q = 0, one gets then 

_ Wit „ , 9 

<?2 = K——Mc. 

G 

4 Conserved mixed currents 

As discussed in the Introduction, and remarked at the end of the previous section, if there 
are non-gravitational fields present in the spacetime then they may also contribute to the total 
superenergy tensor. One important question is whether or not the total s-e current thus con- 
structed may be divergence-free. Of course, this depends on the field equations for the particular 
matter fields involved. In |291l28j . it was proved in full generality that these divergence- free mixed 
s-e currents can always be constructed for the simple case of a minimally coupled scalar field 
(whenever there is a Killing vector in the spacetime). We now summarize here the main results 
in |29| I28j concerning the massless case, which is the one relevant for our present purposes. 

Let <f> be a massless scalar field. The s-e tensor for the scalar field is the basic s-e tensor 
T a pxfiV^[i] V*[i]0} for W(J) and can be easily constructed following the general definition of [2*H] . 
The result is (see [23113123120]): 

-<MV A VVV M V^ - g^VaWVpVrf + 2^^V CT V P 0V CT V^ (32) 
from where we immediately deduce 
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One can straightforwardly compute the divergence of the scalar field s-e tensor 

V a S a px» = 2V /3 V (A 0i^ - gx^V/sV^V^ - V ff 0(2V V (A 0i?% )fta + <?a mj R%/3t V"V» 

so that we realize that the s-e tensor (|32|) is divergence-free in flat spacetime, that is to say, in 
the absence of gravatational field. This allows to construct conserved currents for the scalar field 
in flat spacetime. They are built in a similar way to those formed with the energy-momentum 
tensor or the Bel currents (J5J). Following |291 1281 13T7] . one can define the scalar- field s-e current 
with respect to any three Killing vector fields £ l5 £ 2 > £3 as 

Nevertheless, these currents are not divergence-free in curved spacetimes. In other words, they 
are not divergence-free if one takes into account the gravitational field created by the scalar 
field. This was to be expected: if the s-e concept is to have any physical meaning at all, then 
the total s-e currents, involving both the gravitational and the scalar fields, are the ones to be 
conserved. And this was actually proved in |291l28j : if the Einstein-Klein-Gordon field equations 
are satisfied, that is 

R„u = V^V^ V M V> = 

then the sum of the Bel current and the scalar- field s-e current, which will be written as 

Mi, 6, fs) = Mi, f 2 , fs) + Mi, f2, fs) 
is a divergence- free vector field: 

V Q J a (f 1 ,f 2 ,f 3 ) = 0. 
Observe that J is explicitly defined by 

JfJ. I25 13) = { B (a/3X)^, + S (af3X)tJ,) Cl^.2^3 ■ 

These are mixed conserved currents, and they lead to the conservation of mixed quantities, 
containing both gravitational and scalar fields contributions. Moreover, note that, in general, 
none of the two single currents 2 nor j are divergence-free separately. 

We now want to produce some explicit examples of these conserved mixed currents. To that 
end, we need an explicit solution of the Einstein-Klein-Gordon equations. There are several of 
these in the literature — note that sometimes they are considered as "stiff fluid" solutions. — 
The problem, however, is that most of the solutions have either hypersurface-orthogonal Killing 
vectors, or spherical symmetry, or a GVgroup acting orthogonally transitively, so that the results 
of section |21 apply. This means that, for these particular solutions, the Bel and scalar- field s- 
e currents are, actually, divergence-free on their own. And thus, they do not give the mixed 
conservation we are seeking. What one needs, therefore, is a solution belonging to class A(ii) 
in Wainwright's classification of G2 solutions [321 131j . that is to say, a solution with a non- 
orthogonally transitive G2 group of motions. 4 

Remarkably and fortunately, there exists some explicitly known solutions for spacetimes 
belonging to class A(ii) and with a minimally coupled massless scalar field <j) as source. These 

4 A spacetime with just one non-hyper surf ace- orthogonal Killing vector will also do, but we do not know of any 
such solutions for the Einstein-Klein-Gordon system. 
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were studied by Wils [31] who, among other metrics of the mentioned sort, found the one- 
parameter family of metrics that we are going to use here. 

In general, the line element for such solutions will take the following form: 

ds 2 = -F$dt 2 + F\dx 2 + F 2 [F$dy 2 + F^ 2 (dz + F 4 dx) 2 ] 

where Fq, F±, F 2 and F3 and F4 are functions of t and x only. Wils found the mentioned exact 
solutions making separability assumptions on the functions involved, and the solution we are 
interested in has 

F (t,x) = Fi(t,x) = sinh 2+A I± e &£±L* , 

V6 

F 2 {t,x) =smh^=e? Lx , 
V6 

„ , , 1+2 A Lt (3+A) . 

F 3 (t, x) = sinh"— — e-Sr^ , 



/ \ Lt / 4 A 2 a, 

F 4 (t,x) = cosh-^= yj2- — e*' 



and 



\/4 - 2 A 2 log ^tanh J . 



where L is a positive constant defining the scale and A is a parameter subject to A 2 < 3/2. The 
range of coordinates is constrained by the singularity at t = 0, and thus we take t > 0. 

Obviously, £ = d/dy is a Killing vector that satisfies £ A d£ = 0, so that it is hypersurface 
orthogonal. From Result 12.11 it follows that the Bel current is parallel to £ and hence 
conserved. In turn, and because the total mixed current J(q) is also conserved, the current 
will be divergence-free too. Thus, the sought s-e interchange cannot be found here. We include 
the expression for J(£) for completeness: 

= 4 (A 3 - A + 3)(2A + 3) sinh 2 ^ - A 4 + 20A 2 + 72A + 6s] smh- e ^^e~^ Lx £ 
72 \ y6 / y6 

However, we can also use the second Killing vector, which is if = d/dz, and one has 

rj A drj = -a/3-2A 2 e^r 1 L:c sinh 5+4A ^LdtAdxA dz, 
3 

meaning that 77 will only be hypersurface orthogonal in the extreme cases A = ±-v/3/2. Never- 
theless, it can be checked that the current j(£,ff), as well as 3(£,ry), is in fact proportional to £, 
so that both of them are divergence-free separately once more. Their sum reads 



1 



Lt 



3 (A 4 + 64A 3 + 140A 2 + 40A - 44 - 2(3 - 2A 2 ) 2 ) sinh 4 

V f) 



-4 (14A 4 - 21A 3 - 97A 2 - 30A + 36) sinh 2 

V6 
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The only left possibilities for a true exchange of s-e properties are to be found in the remaining 
conserved currents J(ff) and J(if,£)- And this is indeed the case. Setting x° = t, x 1 = x, x 2 = y 
and x 3 = z, the expressions for their non-zero components are the following: 

J°(V, = j^(3 + 2A) V3^V (4(1 - 2A 2 ) + 3 (3 - 2A 2 ) sinh 2 i0 x 

smh -(ll+6A) U e -3(l+X)L x 

J\ff,0 =-^=V3-2A 2 (^2(3 + 2A) (l - 2A 2 ) + 3 (1 + A) (3 - 2A 2 ) sinh 2 ^0 x 

cosh^lsinh- 6(2+A) ^i e - 3(1+A) ^ 
L 4 ( 

J 3 (r/ ; = ( 204 - 72A - 612A 2 - 320A 3 + 93A 4 + 64A 5 + 4 (81 + 9A - 181A 2 - 

-89A 3 + 50A 4 + 28A 5 ) sinh 2 + 6 (3 + 2A) (3 - 2A 2 ) 2 sinh 4 x 

y6 v6/ 

, «, _i_\\ Lt (9+8A) 
smh -6(2+A) g 

j0 (rf) = ^ Av ^T^ (l - (3 - 2A 2 ) sinh 2 -^=) sinh- 9 -- i£ e -^±Z^ 



ji(^) = __^7^T^ f 2 - (3 - 2A 2 ) sinh 2 ^ cosh ^ sinh 2 5 ' ; ^ > 
18^6 V v^/ \/6 
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(3+7 A) 

e 3 



L 4 f Lt 

J 3 (ff) = 228 + 216A + 44A 2 -3A 4 + 4 (18 + 9A + 6A 2 + 9 A 3 + 2 A 4 ) sinh 2 — - 

216 V \/6 

-2 (3 - 2 A 2 ) 2 sinh 4 ^=V nh ~ 2(5+A) ^= e -d+ 2 ^. 

\/6/ V6 

It can be checked by an explicit computation that these vector fields are divergence- free. We 
want to stress that the corresponding pure currents j(if) and 2(ff), or j(rf,^) and 2(ff, £), are 
not divergence-free in general (unless, of course, 3 — 2A 2 = in which case {£, 77} generate an 
orthogonal transitive G2, in agreement with our previous results). More importantly, we must 
remark that these currents J (if) and J(if,£.) are non-trivial in the sense that they are not linear 
combinations of the Killing vectors. In other words, they would be very difficult to find if we 
did not know about the superenergy concept. The explicit expressions for superenergy tensors 
of the gravitational and scalar fields are essential here. 

A door obviously open by our work can be stated in the form of the following yet unanswered 
question: are there similar non-trivial mixed superenergy currents for general Einstein-Maxwell 
systems? 
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A Appendix 

The 2-form associated with the two Killing vectors spanning the orbits will be denoted by 

w = £ A rj. (33) 



We also define the one-form r y£,A,£,B,£,c)i relative to the vector fields £,a,£,b,£,c, whose com- 
ponents are given by 

r p Cb, fc) = & t% Ra p X a fc R P P ^ (34) 

where we have introduced the indices A,B,C = 1,2 to denote the two independent Killing 
vectors generating the Gi group: £i = £ and £2 = V- ft is convenient to define t (J;a,£b,£cj a s 

V (U, (B,fc) = ^Vau P1P ^Bf3V PUaia ^bRx aP1P2 Ra 1 a 2 a l „ (35) 



where rj a p^ u is the volume element, so that using @, (|3*1)) and ()35[). a straightforward calculation 
allows us to get 

3 (&£a) = \ l-r (Cb,{a,U) ~ 2r (u,U,(b) - ~t (u,Ca,Cb) - \t (Ca,£b,£a 

--j^A R a(3\u ij^Rp^ €b~ Z,B R a(3\a ^Rp^ %A + 

+ l(£ A ■ U)R pPX °R P px„ €b + | (6 • iB)R p/3Xa R Pl 3x, £ A } • 

Clearly the exterior product j Aw contains only terms of the form r Aw and t Aw. We define 
now the 1-forms Jja an d ft a by 

n A ^ = *W a ^ARaf3X p . 

Given that £4 are Killing vectors, we have 

i°i } Rv pa( i = V a [ff,g\p + 2V p r ?[a V^ Q] . 
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Contracting this identity with a third Killing vector, and making the exterior product with w, 
it is straighforward to show 

£U[ a uV] = ^A^VpW^a-p] -^A VP ^ w [^V a ;p] +Ca V /3[ 7 ?^1[q u 'H + ^iA-[v,i]<iw aflu (36) 
Using © it is easy to show that 

^L[a-"V] = 3i ( d * (™ A d Cl) A ^ - *{w A d£ A ) du>) a ^ 

—2 * (?7 A d^ A )^ [pa ^. M] + 2 * (£ A d^)^^^ . (37) 

Bearing in mind these identities, the following definitions will be useful later: 

t Aa = -^a^w^w^ = -^j S A<7 * w ae * w ae , (38) 
&Aa = -O.AlaW^W^ = ~—^Aa * W a€ * W ae . (39) 

In order to obtain an expression for r A w in terms of the above quantities we will use the 
decomposition of the identity 5^ = P®p + Pn~p into its part on the orbit and the corresponding 
orthogonal element at each point, i.e. 

Furthermore, it is convenient to compute first *(r A w) and apply the Hodge dual to the final 
identity. One can then write 

*(r (&, (b,{ c ) a w y = & e B R/x a (Pi P + P\\ P ){Pia + if * ™ ai - 

Taking into account that P^_ p z €Ct * ro" 7 = ^P]_ p z ea * w ae for any 2-form z, this identity can be 
arranged onto its dual form as follows 

r[ a (£a,€b,€c) V] = ^2 {-6RB P A a &Cp * W a[a W pu] 

~ (€b ■ £c) * W Xu Ri2\a(*^A)apu ~ 6 ^A^W^tssRc 5 e[ a W pu ] 

~ (£b ■ £c)Rl212 & (f7A S l[a«V] - CA S 2[a^H)} ' ( 41 ) 

where we have used the indices A, B,C, 1,2 in some places also to denote contractions with 
Ca, £b> £ci £) V respectively, in an obvious manner. 

Let us proceed now with t Aw. Using the previous decomposition of 5^ for the summation 
index v in (|35|) and arranging terms conveniently, we get 

t[a (€i,CB,&) *V] = fF {-^iAWp 2 \iB(3W a2l5 Rc aplP2 Rp ia2 a[ a W^] 

+(U ■ iB)w p ^&R PlP2 x°w a ^R aia2r7[a w pu] 
HU • £b) * w^icRp^x* * w a ^R aia2a[a w pv] ) . 



18 



The first term has the form ry a \£c-, v~a, vhj w^v] with va p = £,A\w px , and hence can be expanded 
using (|41|) . For the second and third terms we still have to apply the same procedure as above 
to the summation index a. After some calculation we can obtain the following expression 



Ha(V\^1p ~ £\^2p)£B~/W T "<£c6W t5 Rr P e[ocW^ u ] + 

+^(U • is) {3\(2t c T w x ° + n c T * w x ")R XaT[a w H 
+ [4(*t c ) af ,u - &(VeXi [a w M - &E 2[a t0 H )] * w Xa R 12X a} ■ (42) 
Using (|41|) and 1)421) we are now ready to obtain the expressions for 



3 {CbXa) A w =~\ ' 



r [€b,U,U) Aw + 2r [£a,U,€b) A io 
+\ t (u, U, A w + - 1 (u, (b,(a) A w 

which, after using £,A*yW Xl (vx^ip ~ £\&2p) = — W&Ap and analogous identities, can be explicitly 
written as (|14jl. 
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